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Abstract. In practical purposes for some geometrical problems in com- 

j^ . puter science we have as information the coordinates of some finite points 

in surface instead of the whole body of a surface. The problem arised here 

is:" How to define a distance function in a finite space?" as we will show the 

OQ | appropriate function for this purpose is not a metric function. Here we 

try to define this distance function in order to apply it in further proposes, 

specially in the field setting of transportation theory and vehicle routing 

problem. More precisely in this paper we consider VRP problem for two 

tyj ■ dimensional manifolds in R 3 . 

O 



oo 



X 



1. Introduction and Preliminaries 
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en ! There are many algorithms for solving problems such as TSP, VRP. All these 

algorithms solve these problems in the case that surfaces are flat Euclidean, 
but in general these problems must be proposed in two dimensional manifolds. 
Here we try to consider the generalized situation. In practical cases there 
are finite points in three dimensional space R 3 . The strategy is to consider 
a two dimensional surface based on these points, Later we define function 
for the distance between these points in the defined surface. To do this first 
we define a convenient two dimensional manifold respect to these points, and 
in the second step we define the distance between two points respect to this 
manifold. 

For the first level least square method have been generalized to define con- 
venient two-dimensional manifold. For the second level,it is expected to define 
distance between two points by finding vertical feet and then geodesic which 
connect these vertical feet. Although on the surface this definition seems a 
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natural definition, there are examples which show this definition is not an ap- 
propriate one. Therefore we try to define a new suitable definition. In the new 
definition for every two points A, B we consider a two- dimensional manifold 
which passes through given two points and also it satisfies the least square 
method for the rest of points, by applying this we define the distance between 
A, B. We consider the length of the geodesic which is passed through two 
points A, B on the derived surface. Further, we show this definition is better 
than the former one. In the first algorithm by considering the natural number 
as fixed, we passed through these points a polynomial(manifold) degree n, like 

n—i n 

i (n—d—i) 



'1-1) z = ^2^2 a i,(n-d-i)X t y 



d=0 i=0 

In the second algorithm to find the distance between two points A, B as 
we have seen in the first algorithm, The manifold should met the following 
requirements: 

(1) The manifold passes through these two points. 

(2) This manifold is chosen such that it satisfies the least square conditions for 
the rest of points. 

Let we consider the following example in order to show that why we prefer the 
second definition. 

Example 1.1. Consider a finite set of points {Pi, P%, ..., P n } in R 3 such that 
{Pi, P2, ..., P n -i} is encompassed by a circle C(0, r) in the xy -plane and P n is 
P n = (0, 0, d) with d ^> r. Regarding the first definition distance between Pi 
and P n is less than r(for large enough n ), but in accordance with the second 
definition distance between Pi and P n is close to d as we expect intuitively. 

The second definition is not a definition of metric as we see later. In the 
next section,we describes process to fined the distance between each pair of 
points in the finite set of points. 

2. GENERALIZED LEAST SQUARE METHOD 

There exist numerical methods for obtaining a useful two dimensional man- 
ifold which can be appropriate for fitting points Pi, P2, ..., P& in M 3 . We would 
like to find a surface of degree n can be fitted by Pi, P 2 , ..., P^.At the first step 
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we write complete form of equation of degree n as follows: 



n—i n 

EE 

d=0 i=0 



"i,(n—d—i)X V 



i (n—d—i) 



At the second step we consider the sum of square differences of this equation 
with the given points , namely: 
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At the third step we find parameters such that this sum would be mini- 
mized. Then by derivating with respect to these parameters, and equaling it 
to zero we obtain the following equations. 



{ 
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for i — 0, .., n; d — 0, ...n — i.} 



By solving this system of equations we have a polynomial manifold. This is 
an inherent generalization of least square method to find a line passing a set 
of points. It is worth mentioning to note that in solving this system there is a 
matrix with the possibility of having the determinant equal to zero,in the case 
that the determinant of matrix is equal to zero we have infinite surfaces that 
can be fitted with P 1; P 2 , ..., P^. 

Example 2.1. For points Pi,P 2 , ...,P±o such that P$ = (xi,yi,Zi) , finding a 
surface using generalized least square method leads to a matrix with determi- 
nant equal to zero. Therefore we have infinite classes of parameters to define 
a surface. 



Table 1. The Coordinate of Points 
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3. Distance function 

Definition 3.1. Let X be a finite set of points in IR 3 , for any nGN, distance 
function d n : X x X — > IR has been defined as follows: 

We consider x,y G X ,d n (x,y) is the length of geodesic which has passed 
through {x, y} on the surface of degree n , and satisfies the condition of least 
square for the rest of points. 

At first it must be showed that this function is well defined. i.e, for any two 
points x, y there is a unique surface which passes through points and satisfies 
least square condition for the other points, Uniqueness of the surface is related 
to the determinant of the equation of the system for least square, i.e, if deter- 
minant is non-zero then related system has a unique solution for coefficients 
of surface, therefore we have a unique surface with desired conditions. But if 
determinant is zero, there are infinite surface that satisfy conditions. We show 
in the next theorems we are able to define the distance function as desired. 

According to this definition, 2 dimensional manifold of n degree which has 
passed two points Pi, P2 should be obtained at first (we want to obtain length 
of this points.) and with respect to other points satisfies the least square con- 
ditions. The type of equations will change in comparisons with the mentioned 
method. In the new equations, two primary equations show fitness of Pi,p2 
on the surface,the remaining equations should have Maximum fitness in the 
respect to other points. Thus we have the following equations : 

En—i \r^n j (n— d— i) 
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Now we compute the length of geodesic passing through two points (Pi,P2). 
Let Mc(Pi, ..., Pfc) be matrix associated to system equation Q dD — = 0. One 
possible problem during solving the equation system is the possibility of zero 
determinant of Mc(Pi, ■■■, Pt 



k 



Theorem 3.2. Consider a set of points X = {P , Pi, ..., P^}, assume M c (Pq, Pi, 
P2, ..., P n ) was equation system of least square method. For any Pq G X there 
is sequence of points P 1; P 02 , P 03 , ... such that linij^ooPQ^ = P and for any 
k G N , determinant ofY^ = {P ok , Pi, P2, ...,P„} was non-zero. 
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Theorem 3.3. If C is a surface such that limj_ s . 0O Ci = C , then the limit of 
geodesies g^ which passed through between to points P n , P m is g,i.e Hindoo ^ = 

9 

Remark 3.4. d have been obtained by least square methods and Ci are poly- 
nomial surfaces ( 11. ip . condition lim^oo Ci = C is equivalent to lim^oolafcjj = 

Ck,l 

Let C be a manifold and the determinant of related matrix to C has zero 
determinant i.e det(Mc) = 0. Now we consider {Cj}?^ such that all of these 
surfaces passed through x, y and the determinant of the associated matrix of 
them is not zero. They define the sequence of geodesies {<?i}^i- By previous 
Theorems lim^oo^ = g such that g is the geodesies of C and lim^oo ||^j|| = 
\\g\\. So distance function is well defined. 

Definition 3.5. We are capable to generaliz the previous definition to dx '■ 
M. 3 x M. 3 — y K such as for any x,y G IR 3 ,we define d as follows: dx(%,y) = 
dxu{x, y }(x,y). 

we call (M 3 ,<i) 3-finite distance space. In a natural way we are able to 
generalized the above theorem to M. n . we call (M ra ,<i) n-finite distance space. 
In this paper we succeed to define the distance function to find distance be- 
tween points in two dimensional manifolds. Later on, we will apply it to some 
problems in transportation theory. 

Lemma 3.6. n-finite distance spaces are not metric spaces. 

The following counter example shows that d is not a metric. 

Example 3.7. Let {Pi, ..., Pio} be points as show in Table l.For points Pi, P3, P10 
d(10, 1) = 18.84342860, d(10, 3) = 12.78704132, d{l, 3) = 4.123340349 
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